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Abstract 



Higher order digital nets have recently been recognized as one of the 
most promising branches of quasi-Monte Carlo methods. The notable fea- 
ture of higher order digital nets is that they can exploit the smoothness of 
a function for numerical integration and achieve the improved convergence 
rate of the integration error for smooth functions. One prominent con- 
struction of such nets is based on digitally interlacing the components of 
the classical digital net whose number of components is a multiple of the 
dimension. In this study we consider the weighted unanchored Sobolev 
spaces of high smoothness and derive an upper bound on the mean square 
worst-case error for digitally shifted higher order digital nets. Employing 
our obtained bound as a quality criterion, we prove that the component- 
by-component construction can be made efficient use of to obtain good 
polynomial lattice rules that are used for interlaced components. Through 
this approach we are able to get some tractability results under certain 
conditions on the weights. Numerical experiments confirm that the per- 
formance of our constructed point sets often outperforms in terms of our 
introduced quality criterion the performances of higher order digital nets 
with SoboP sequences and Niederreiter-Xing sequences used as interlaced 
components, indicating the usefulness of our algorithm. 

1 Introduction 

In this paper we study the multivariate integration of smooth functions defined 
over the s-dimcnsional unit cube [0, l) s , 



Quasi-Monte Carlo (QMC) methods approxiamte /(/) by averaging function 
evaluations at uniformly distributed TV points xq, . . . , jcjv— i, 
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We refer to [10, 18] for general information on QMC methods. The Koksma- 
Hlawka inequality ensures that QMC rules typically achieve the worst-case error 
of order N~ 1+s for any <5 > when the integrand has bounded variation in the 
sense of Hardy and Krause. Two prominent ways to generate a uniformly dis- 
tributed point set are lattice rules, see for example [20, 27], and digital (t, to, s)- 
ncts, see for example [13, 20]. It has long been known that it is possible to 
achieve the improved convergence rate for periodic smooth functionss by using 
lattice rules, whereas neither rule can exploit the smoothness of a non-periodic 
function so as to achieve the improved convergence rate. 

Recently, higher order digital nets have been proposed in [7, 8] and shown to 
achieve the improved convergence rate for numerical integration of both periodic 
and non-periodic smooth functions. Numerical integration using higher order 
digital nets in a reproducing kernel Hibert space consisting of smooth periodic 
functions was studied in [7] , and that in a normed Walsh space consisting of all 
smooth non-periodic functions was analyzed in [8]. Later it was shown in [2] 
that the result from [8] can be achieved in a reproducing kernel Hilbert space 
setting that is different from what is studied in [7] . More specifically, a weighted 
unanchored Sobolev space of smoothness a G N, where N is a set of positive 
integers, was considered as a reproducing kernel Hilbert space and it was proved 
that higher order digital nets achieve the optimal convergence rate not only of 
the worst-case error but also of the mean square worst-case error with respect 
to random digital shift. We note that the case with a = 1 has been extensively 
studied in the literature already, see for example [13, Chapter 12]. More recently, 
higher order scrambling was studied in [9, 15] to achieve the (nearly) optimal 
rate of the root mean square error for smooth functions. 

There are two construction algorithms of higher order digital nets so far as 
known. One has been studied in several papers such as [3, 4, 12], which makes 
use of polynomial lattice rules introduced in [21] by generalizing the original 
definition. These rules are referred to as higher order polynomial lattice rules. 
Regarding explicit constructions of higher order polynomial lattice rules which 
achieve the nearly optimal rate of convergence, the component-by-component 
(CBC) construction has been studied in [3, 4]. Even with efficient calculation 
of the worst-case error implemented in [4], however, a computational cost of 
0(sdN d log TV) operations using 0(N d ) memory is required for d G N. Whereas 
higher order polynomial lattice rules thus constructed can achieve the worst-case 
error of order TV" min (<M)+<5 for functions in a normed Walsh space of smoothness 
a, the exponential dependence on d of the number of operations degrades the 
availability of these rules as d increases. 

The other algorithm for constructing higher order digital nets has been pro- 
posed in [7, 8], which applies a digit interlacing function to the components of 
the classical digital net whose number of compenents is a multiple of the di- 
mension. Hence the major advantage of this algorithm is that one can make 
use of the existing digital (t, m, s)-nets including those of Sobol' [29], Faurc [14], 
Niederreiter [19] and Niederreiter-Xing [22]. Suppose we have digital (i', to, ds)- 
nct in base b for d € N. We denote by y n = (y n ,i, ■ • ■ , Un,ds) the n-th point of 
this net and by y n j = y n j,i^ _1 +Un,j,2b~ 2 + ■ ■ ■ the 6-adic expansion of y n j for 
1 < j < ds. Then the n-th point x n = (x n i, . . . , x r ^ s ) S [0, l) s of higher order 
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digital nets consisting of b m points is given by 

oo d 

X n ,j = fd(2/n,d(j-l)+lJ ' • ■ ^Vn,dj) ■= y n A(j-l)+r,ab~ r ~^ a ~ 1 "' d , 

a—1 r—1 

for 1 < j < s. Here we refer to T>d as a digit interlacing function of factor 
d. The disadvantage of this construction algorithm is that the i-valuc of dig- 
ital (t,m, s)-nets increases at least linearly in s, so that it becomes hard to 
obtain a bound on the worst-case error independent of the dimension. This 
observation motivates us to replace the existing digital (t',m, ds) -nets by suit- 
ably chosen polynomial lattice rules in dimension ds that are used for interlaced 
components. In the similar context, there exists a successful result in [15] where 
scrambled polynomial lattice rules are used as interlaced components to con- 
struct higher order scrambled digital nets. It was shown that we are able to 
obtain a better dependence on the dimension of the root mean square error. 
Moreover, as compared to deterministic higher order polynomial lattice rules, 
the computational cost for the CBC construction could be significantly reduced 
to 0(sdN log N) operations using O(N) memory. Thus, it is quite natural to 
think of using polynomial lattice rules for interlaced components to achieve the 
nearly optimal rate either of the worst-case error or the mean square worst-case 
error with respect to some randomization while obtaining a better dependence 
on the dimension. 

In this study, we consider the weighted unanchored Sobolev spaces of high 
smoothness as studied in [2] . We derive a computable upper bound of the mean 
square worst-case error in this space setting for digitally shifted higher order 
digital nets. Employing this upper bound as a quality criterion, we prove that 
the CBC construction can be made efficient use of to obtain good polynomial 
lattice rules that are used for interlaced components. In the next section, we 
introduce the necessary background and notation including Walsh functions, 
polynomial lattice rules, higher order digital nets and their randomization, and 
the weighted unanchored Sobolev spaces of high smoothness a S N. In Section 
3, we study the mean square worst-case error in this space setting for digitally 
shifted higher order digital nets with an aim to introduce a computable upper 
bound on the error. We show in Section 4 that the CBC construction is applica- 
ble to obtain good polynomial lattice rules as interlaced components such that 
digitally shifted higher order digital nets thus constructed achieve the nearly 
optimal rate of convergence and that we are able to get some tractability re- 
sults under certain conditions on the weights. We remark that the fast CBC 
construction using the fast Fourier transform as introduced in [23, 24] is avail- 
able in our current setting. This enables us to proceed the construction with 
0(sdN log N) operations using O(N) memory. Finally, we conclude this paper 
with numerical experiments in Section 5. 

2 Preliminaries 

Throughout this paper, we use the following notation. Let N be a set of positive 
integers and No := NU {0}. Given a prime b, let F& := {0, . . . , b — 1} be the 
finite field consisting of b elements. The operators © and denote the digitwise 
addition and subtraction modulo b, respectively. That is, for x,x' G [0, 1) with 
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6-adic expansions x = YliLi 1 an< ^ x ' = x 'i^ 1 wnere x ii x '% £ IF& , © 

and are denned as 



where yi = Xi + x[ mod 6 and y[ = Xi — x\ mod 6, respectively. Similarly, we 
define digitwise addition and subtraction for non-negative integers based on 
those 6-adic expansions. In case of vectors in [0, l) s or Ng, the operators © and 
G are carried out componentwise. We call x £ [0, 1) a 6-adic rational if it can 
be represented by a finite 6-adic expansion. Further we shall use the notation 
:= {1, . . . , a} for a e N for simplicity. 

2.1 Walsh functions 

Walsh functions were first introduced in [30] for the case with base 2 and were 
generalized later, see for example [5]. We refer to [13, Appendix A] for gen- 
eral information on Walsh functions. We first give the definition for the one- 
dimensional case. 

Definition 1 Let b > 2 be an integer and ujb = e 27ri / fc . Let us denote the b-adic 
expansion of k £ No by k = kq + K\b + • • • + K a -\b a with Ki <G {0, . . . , b — 1}. 
Then the k-th b-adic Walsh function bwalfe : [0, 1) — > {1, Ut,, . . . , w^" 1 } is defined 



for x £ [0, 1) with its b-adic expansion x = x\b 1 + x^b 2 + • • • , that is unique 
in the sense that infinitely many of the Xi are different from 6—1. 

This definition can be generalized to higher-dimensional case. 

Definition 2 For dimension s > 2, let x = (x±, . . . ,x s ) € [0, l) s and k = 
(k\, . . . ,k s ) € Ng. Then the k-th b-adic Walsh function fcwalfc : [0, l) s — > 
{1, u>b, ■ ■ ■ , * s defined as 



Since we shall always use Walsh functions in a fixed base 6, we omit the subscript 
and simply write walfe or wal& in the remaining of this paper. 

2.2 Polynomial lattice rules 

Given a prime 6, let us denote by F&((x _1 )) the field of formal Laurent series 
over Fft. Every element of Fb((x -1 )) has a form 




as 



b wal k (x) 



Xl^O H \-X a K a -i 

'b 



s 



b w&\ k (x) = 6 wal fcj (xj). 

3 = 1 



l = t\x 



l—w 
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where w is an arbitrary integer and all ti € Ft>. Further, we denote by ¥i,[x] the 
set of all polynomials over F&. For a given to € N, we define the mapping v m 
from Fh((a;~ 1 )) to the interval [0, 1) by 

(OO \ TCI 

l—w / I— max(l,u>) 

We often identify a non-negative integer k whose o-adic expansion is given by 
k = KQ + Kib+- • • + K a _i6 a_1 with the polynomial over ¥i,[x] as k(x) = ^0 + ^1^+ 
• • • + Ka^x"- 1 . For k = (k u ...,k 8 ) £ (¥ b [x}) s and q = ( qi ,...,q s ) £ (¥ b [x]) s , 
we define the inner product as 

s 

k-q E k jH e V b[x], 
i=i 

and we write q = modp if p divides q in F&[:r]. Using these notations, a 
polynomial lattice point set is constructed as follows. 

Definition 3 Let to, s S N. Lef p € F&fx] 6e an irreducible polynomial with 
deg(p) = to and let q = (qi, . . . , q s ) £ (Fh[x]) s smc/i t/iat ^ and deg(qj) < m 
/or 1 < j < s. A polynomial lattice point set is a point set consisting ofb m points 
such that 

(^))^>- 

for < n < b m . A QMC rule using this point set is called a polynomial lattice 
rule with generating vector q and modulus p. 

In the remaining of this paper, we denote by Qb m ,s{q,p) a polynomial lattice 
point set, implicitly meaning that deg(p) = m and the number of components 
in a vector q is s. The concept of the so-called dual net plays an important role 
in the subsequent analysis. 

Definition 4 Fork G N with its b-adic expansion k = ko + ki& + - • • + K a _io a_1 , 
we denote the truncation of the associated polynomial k(x) S ¥{,[x] by 

tr m (fc) = K + K lX + ■ ■ ■ + Kn^a,™)-!!"'"'" 1 '" 1 . 

For a polynomial lattice point set Qb m ,s(Q,p), the dual net is defined as 

D x {q,p) = {k = (h, ...,k s )€N s : tr m (fc) • q = modp}, 

where tr m (fc) = (tr r „(fci), . . . , tr m (k s )). 

Futhermore, we shall use the following two lemmas in this paper. The first 
lemma bridges between the dual net D ± (q,p) and Walsh functions. The proof 
is straightforward from the above definition of D J -(q,p), [13, Lemma 10.6] and 
[13, Lemma 4.75]. The second lemma states about the group structure of a 
polynomial lattice point set Qb™,s(q,p)- The proof is also straightforward from 
[13, Lemma 4.72]. 



n(x)qi(x) 
p(x) 
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Lemma 1 Let Qb m ,s(q,p) be a polynomial lattice point set and D ± (q,p) be its 
dual net. Then we have 

J-E 1 wal fc(a;n ) = (i **^(«.*). 
um A^i k\ n> 1 n otherwise. 



n=0 

Lemma 2 For a prime b, any polynomial lattice point set Qb"\s(9,p) is a sub- 
group o/(F™ s ,©). 

2.3 Higher order digital nets 

QMC rules based on higher order digital nets can exploit the smoothness of an 
integrand so that they achieve the optimal convergence rate of the integration 
error for functions with smoothness a £ N. This result is based on a bound 
on the decay of the Walsh coefficients of smooth functions [2, 8]. Explicit con- 
struction of higher order digital nets by means of a digit interlacing function 
was given in [7, 8]. Although we have already mentioned about that function 
in the previous section, we describe the interlacing algorithm in more detail in 
the following. 

Since the interlacing is applied to each point separately, we use just one 
point to describe the procedure. Let y £ [0, l) ds with y = (yi, ■ ■ ■ ,Uds) an d 
denote the 6-adic expansion of each coordinate by yj = yj^b^ 1 + yj^b^ 2 + • ■ • , 
that is unique in the sense that infinitely many digits are different from 6—1. 
We then obtain a point x £ [0, l) s by interlacing the digits of non-overlapping 
d components of y in the following way: Let x = (x±, . . . , x s ) where 



-(a-l)o 



for 1 < j < s. We denote this mapping by T> d : [0, l) d — > [0,1) and we 
simply write xj = T> d (y d{j _ 1)+1 , . . . , y dj ). Further we write x = V d (y) := 
{V d (yi, ...,y d ),..., Z>d(j/d( s -i)+i, • • • ) Vds)) when x is obtained by interlacing 
the components of y. Thus it is obvious that in order to construct a good 
higher order digital net consisting of N points in [0, l) s we need to carefully 
choose N points in [0, l) ds . 

In this paper, we are interested in applying polynomial lattice rules to gen- 
erate a point set in [0, l) ds that serves as interlaced components for higher order 
digital nets. For clarity, we give the definition of higher order digital nets based 
on polynomial lattice rules. 

Definition 5 Let b be prime and m,s,d £ N. Let p £ F&[a;] be an irreducible 
polynomial with deg(p) = m and let q = (qi,...,q ds ) £ (¥i,[x}) ds such that 
qj =/= and deg(qj) < m for 1 < j < ds. We construct a higher order digital net 
consisting of b m points {xq, . . . , x^m^i} C [0, l) s as follows. The n-th point x n 
is obtained by 

x n = V d {y n ), 

where the point y n £ [0, l) ds is the n-th point of a polynomial lattice point set 
Qb m .ds(Q-,p) which is given as 

( (n(x)qx(x)\ (n{x)q ds [x) 



p{x) J ' ' V P( x ) 



() 



for < n< b m . 

In this construction algorithm, the search for good b m points in [0, l) rfs to be 
interlaced has now been reduced to finding good polynomials q = (q±, . . . , qda), 
which is the particular interest of this paper. 

Randomization of QMC point sets is useful to obtain some statistical es- 
timate on the integration error. Especially for randomization of higher order 
digital nets, two algorithms have been discussed in the literature. One is digital 
shift, see [2, 7], and the other is higher order scrambling that is a generalization 
of Owen's scrambling introduced in [25], see [9, 15]. Since we are concerned 
with the former in this paper, we follow [7, Section 6] to introduce some basic 
concepts of digital shift here. 

Let Qfem jS = {x ,...,Xbm-i} C [0, l) s with x n = (x n ,i, . . . , x„ tS ). We 
denote the 6-adic expansion of x„ t j by x n j = x n) j t ib + Xnj^b + • • • for 
1 < j < s. Also, let er = (ai, ...,cr s ) where Oj is independently and uniformly 
distributed in [0, 1) for 1 < j < s. We also denote the 6-adic expansion of <7j by 
tjj = Oj^b^ 1 + Gj^b + • ■ ■ for 1 < j < s. Then the randomly digitally shifted 
point set Qb™,s,o- = {zo, z bm _ x } is given by 



where 



Z-n © & \Zn,l ; * • • : ^n..s): 



for 1 < j < s. In order to compute the mean square worst-case error, the next 
lemma is required. We refer to [13, Lemma 16.37] for the proof. 

Lemma 3 Let xi,X2 S [0,1) and let z\,Z2 S [0,1) be the points obtained by 
applying an i.i.d. random digital shift to X\ and X2, respectively. Then we have 

walfe(xi 9^2) if k = I, 
otherwise. 



E[wal fe (zi)wal/(z2)] 



2.4 Weighted Sobolev spaces of high smoothness 

Here we follow [2] and introduce the weighted reproducing kernel Hilbcrt space 
considered in this paper. The concept of weighted spaces was introduced in [28], 
which provides an essential insight into tractability properties of multivariate 
problems from the viewpoint of information-based complexity. Let us start with 
the one-dimensional unweighted case. The inner product of the Hilbert space is 
defined for a £ N by 

</,S>iW> =E f f {r) (x)dx f g^{x)dx + f f^(x)g^(x)dx, 
r=0 Jo Jo Jo 

where we denote by f^ the r-th derivative of / and set /^°^ = /. Let 
ll/llffi,c(i) = y/(/> /)ffi,c(i) be tn e norm of / associated with 
now define the function 

^ B r (x)B r (y) ^ a+1 B 2a {\x-y\) 
A ^W( I 'W = L (H)2 + (2a)\ ' 



We 
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where B r denotes the Bernoulli polynomial of degree r. We refer to [1, Chap- 
ter 23] for information on Bernoulli polynomials. The reproducing kernel for this 
one-dimensional case is given by 1 + K l a ^(x, y). That is, for any / 6 i?i ja m, 
we have 

f(x) = (f,l + #!,„,(!) (■,£))#!,„,(!, • 

We now consider the higher-dimensional weighted case. The inner product 
for the s-dimensional weighted unanchored Sobolev space H s an of smoothness 
aeN and non-negative weights 7 = (j u ) u ci s is defined by 

(f,9)Hs, a ,-y 

= E^ 1 E E 



uCI s 



X 



/ / /C r «\«'""°> (*)<&-„ [ g<- r *» a <"V(x)dx- v dx v , 

./[O,!) 1 " 1 J[0,1) S - 11,1 7[0,l) = -l"l 



where we have used the following notation. For r u \ v = (rj)j£ U \v, we denote by 
(r u \ v , a v , 0) the vector in which the j-th component is Tj for j £ u\v, a for j e u, 
and for I s \u. For u C I s such that j u = 0, we assume that the corresponding 
inner double sum equals and we set 0/0 = 0. Let = y/Jf, /)h s . q , t 

be again the norm of / associated with 7J S:Q/7 . 
The reproducing kernel for H s ^ an becomes 

\- TT (sr B r{xj)B r ( yj ) ^ a+ i B 2a (\xj -Vj\) \ 
~ 2^ T«H I 2^ ( r! )2 +^ L > r 2a )! j ' 

uCI s jeu \r=l v ; K ' J 

where we set 

J'€0 

That is, we have for any / g -ff s , Q .-y 

3 Mean square worst-case error 

In this section, we derive an upper bound on the mean square worst-case error 
in the space H s ^ a ^ for digitally shifted higher order digital nets. To begin with, 
the worst-case error in the space H s ^ an is defined by 

e(Q bm>s ,H Sian )= sup \I(f)-i(f)\, 

ll/llH.,a, T <l 

where an arbitrary point set Qt m , s consisting of b m points is used for quadrature 
points. The initial error is given as 

e(Qo,a,H a ,a,-f) = SUp |/(/)|- 

II/IIh. „,<i 
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From [2, Theorem 13], we have the following simple expressions for the squared 
worst-case error and the squared initial error in the space H s>an 

e 2 {Qb^. Sl H s . an ) = / / K s>an {x,y)dxy 

J[0S)° J[0,1)° 

2 bm " 1 - ' b '"- 1 



6 " "o •'[0,1) 
= _ T0 + ^ X K s,a,-r(x n ,x n ,), 

n,n'—0 

and 

e 2 (Qo,s, H at a t ~f) = 70, 

respectively. Let us now consider a randomization of Qb m , s by using digital 
shift. The mean square worst-case error for digitally shifted point set Q&m Sjer , 
denoted by e 2 (Qb™, s , Hs,a,-y), is defined and calculated as 

^2 



e (Qb m , s > Hs,a,-t) : =E ^(Qt^^^j-ffs^^)] 

j 6»-l 

= ~ 70 + 52^" X E I^.oi.T^n © o", aj„' © er)] 

n,n'=0 

=52^ zj 5i 7«n E [ if iA(i)( i ' i J 9<7 J' ;,: «'j ffl,7 j)] ■ 

n,n'=0 0/«CJ a j£u 

On the other hand, the mean square initial error is still given as 

e 2 (Qo,s,-ff s ,a, 7 ) = 70- 
We now consider the Walsh expansion of K l a ^(x,y). We write 

00 

#i,a,(i)(s, y) = 51 ^i,a,(i)( fc )0 wal fc( a: ) wal Ky)! 

where the (fc, Z)-th Walsh coefficient K la ^(k,l) is defined by 



^i,a,(i)(M) = / / A'i.Q,(i)( x 'y) wal '=( a; ) wal '(y) da;d y- 

Jo Jo 

Using the above Walsh expansion of K-i a !\\{x, y) and Lemma 3, we further 
have on the mean square worst-case error 

^ 6 171 — 1 00 
n,n'=0 0^uCJ s i£ufe=l 

1 b '"- 1 

= zJ 7u X] X] wal (k u ,0)(x n OXn>), (1) 

0/uC7., fc u GNl"l n,n'=0 

where we write K la ^{k) = K lj0l> ^(k, k) for short and -Ki, a ,(i)(fcu) = ]J jeu K la ^(kj). 
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3.1 A bound on Walsh coefficients 

From [2, (13) and Proposition 20], we have a bound on Walsh coefficients 



•^l,a,(l)( fc ) 



where fi a (k) is the weight first introduced in [8], which is defined as 

fj, a (k) := Ol + h a m in(iy,a), (2) 

for k £ N with its 6-adic expansion k = Kib ai ~ 1 + • • • + n l/ b a "^ 1 such that 
< Ki, ■ ■ ■ , k v < 6 and a\ > ■ ■ ■ > a u > 0. We set /i a (0) = 0. Moreover, D a ^ is 
given by 

D a>b = max (C abv + C 2a ^ a -^) , 

l</i<a \ ' / 

in which C f a b v and Cioi.b &re respectively given by 

a 

where C hb = (2 sin( 7 r/6))- 1 and C T , b = (1 + 1/6+ 1/(6(6+ l))) T - 2 /(2 sin(7r/6)) T 
for r > 2, and 



C2a, 



2(1 + 1/6+1/(6(6+1))) 



2a-2 



(2sin(7r/6)) 2Q 

Furthermore, we write [i a (k u ) = J2j^u Ha(kj). We then have 



Inserting this bound on the Walsh coefficient into (1), we have thus far 



b m -l 



e 2 (Q b m !S ,H s , an )< 1» D a!b E b ~ 2 ' la(ku) -g^ E wal (feu ,o)(x„ea ; „0 



n,n'— 



3.2 A bound on mean square worst-case error 

We continue the study of the mean square worst-case error in the space if SiQ ..y to 
derive a computable upper bound. Following [9], we first define a digit interlac- 
ing function of factor d for non-negative integers. For d £ N and k\ , . . . , kd £ Nq, 
we denote the 6-adic expansion of kj by kj = Kj > o+Kj ) ib+- ■ ■ for 1 < j < d. Then 
a digit interlacing function £ d gives one non- negative integer from k\ , . . . , kd by 

d oo 

£ d (k u • • • , M = E E K( j -l)d+r,ab ad+r - 1 - (3) 
r=l a=0 

We extend this function to vectors as 

£d{h, ■ ■ ■ , k ds ) = (£d{h,- ■ ■ , k d ), ■ • ■ , £d(k d ( s -i)+i, • • • , fas))- 
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From the definition of Walsh functions, we have 

d 

i=i 

for any fci, . . . , k d G N$ and (y 1} . . . ,y d ) G [0, l) d . 

We now recall that the point set Qb™, s — {xo, ■ ■ ■ >a?6 m -i} is constructed 
according to Definition 5. Using the above property on the digit interlacing 
functions T> d and £ d , and the group structure of polynomial lattice rules as 
described in Lemma 2, we have 

< E 7,( U )< ( ; }1 E ^^^^"'^E-'M^o^^e^)) 

0/«C7 d3 fc„GNM n,n'=0 

= E ^M^f 5 ' E b-^&^M, (4) 
(fe u ,0)e£>- L (g,p) 

where (fe u , 0) is here a vector in Nq s such that kj G N for j € u and fcj = for 
j ^ u. Furthermore, we have denoted by v(u) the set of 1 < j < s such that 

un{d(7-i) + i,...,<#}^0. 

The next lemma gives a lower bound of (J, a (£ d (k u , 0)), from which we can 
derive a computable upper bound on e 2 (Qb m .ds(q,p), H s>a) y). 

Lemma 4 for $ ^ u <Z I ds and k u G N' u l, we /iai>e 

/i a (^d(fe«,0)) > min(o!,d)y^/ii(fcj) + ^a|«| - ^a(2d- l)\v(u)\. 

jeu 

Proof. We denote by Uj the set of u fl — 1) + 1, ... , dj} for j G v(it). Then 
obviously we have 

A*a(£<i(feu)0)) = ^2 ii a (£ d (k U] ,0) d ), 

j€v{u) 

where (k Uj ,0) d is a vector in Nq such that fc; G N for I G Uj and ki = for 
Z G — 1) + 1, ... , dj} \ uj. Thus it suffices to prove that 

H a {£d{k Uj ,Q)d) > min(a,d) ^ ^i{h) + ^a\uj\ - -a(2d- 1), 

for j G u(u). Considering the definition of the weight fi a (k) as in (2), in order to 
precisely evaluate the weight fJL a ifid{k\, . . . , kd)) in (3) for a given k\, . . . , k d G 
Nq, we need to reorder the summand according to the value of ad + r. Instead 
we give a lower bound by only looking at the most significant digit in the inner 
sum of (3). We denote by {tTj(1), . . . , nj(d)} the permutation of {1, ... , d} for 
a fixed j such that 

(Ml(fcdC7-l)+7r j (l))-l)d + 7Ti(l) > ••■ > (Ml(fcd(j-1)-Hr,(d)) ~ X ) d + TTj(d). 
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Here we note that only the first \uj\ terms are positive since fcafj— i)+7r-(0 = 
for \v,j\ < j < d. If a > \v,j\, it is obvious from the definition of £d{k u ,0) d to 
have 

Kl 

fi a (£ d (k Uj ,0) d ) [(Mi(^d(j-i)+7r 3 (i)) — l) (i + 7Tj(z)] 

i=i 



> min(a, d) ~~ a ^ + — ( | ?^ ^ | + 1). 

Otherwise if a < owing to the averaging argument of the pcrmutated terms, 
we obtain 

a 

H a {£d(k Uj ,0)d) > ^2 [(Ml(fcd(3-l)+7r 3 -(i)) - 1) d + 7T,-(i)] 



^ l77 H [(Mi( fc dO--i)+^-(i)) - 1) <* + Ti(*)] 

1 3 ' t=l 

> a ^ - ad + + 1). 



Putting the two cases above together, the result follows. □ 

Now we arrive at the following theorem on an upper bound on the mean 
square worst-case error for digitally shifted higher order digital nets that are 
constructed using the digitally interlacing of polynomial lattice rules. Since the 
proof is straightforward by applying Lemma 4 to (4), we omit the proof. 

Theorem 1 For higher order digital nets constructed according to Definition 
5, the mean square worst-case error in the space H s>a ^ with respect to digital 
shift is bounded from above by 

e 2 (Q bm<s ,H Sta „)< Mu)D l S )l E r a , d (K), 

(k u ,0)€D ± (q,p) 

where we define D a d = b < - 2d ~ 1 ^ a D Etyd , 

f 1 if* = 0, 

r a ,d{K) ^ & -2mi„( Q ,d) All (fe)-a otherwise) 

and r a4 {k u ) = H j&u r a ^(kj). 

In the remaining of this paper, we employ this upper bound as a quality 
criterion for construction of polynomial lattice rules. For simplicity, we denote 
the bound by 

B a ,d n (q,p) = ^2 lv(u)D^' d )][ ^2 r a,d( k u)- 

(fe»,0)e£'- L (q,p) 



12 



The next corollary provides us with a computable formula of B ai d,y(q,p). 
We write log fc for the logarithm in base b and set &(2min(o ; ,d)-i)|>g i ,0j L 

Corollary 1 Let B a ^,-y(q,p) be given as above. We have 



b m -l 



n=0 B/mC/, 



1 + II C 1 + KAVr>.,d{j-i)M)) 
1=1 



where for y £ [0, 1) 



b _ 1 _ {,(2min(oi,d)-l)Llog i) yj ^2min(a,ii) _ ^ 



^oi^min(a!,rf) _ 5) 

Particularly in case of product weights, that is 7„ = IJjeu 7j' we reduce to 



b m -l 



B a ^( q , P )=-i+± ri 



1 - ljD a>d + 'JjDa^ JJ (1 + 0a,d(j/n,d(j-l)+j)) 

n=0 j = l 

Proof. Because of the property of the dual net D ± (q,p), we have 

1 fem - 1 

B a ,d n (q,p) = 12 7K«)^id E r ",d( fc «)l^ $1 wal (fe«,o)(y„) 



6 m -l 



fe«f 



iE E 7» w j5 )l IIE r MfeM fej (!/ BJ -). (5) 



n=0 B^ M C/ dl j£ukj = l 

We follow the similar lines as in [11, Section 2.2] to obtain for y G [0, 1) 



y"Vq,d(fc)walfc(?/) = 4> a ,d{y)- 



fe=i 



We arrange (5) by collecting the terms associated with a given w C I s such that 
t>(it) = ty. We then have 



6 m -l 



B a4n {q,p) = — E II 



n=0 0^u)C/ s 



1 + II ( X + <?W(2/n,d(j-l)+;)) 



/=1 



Hence the result for the first part follows. For product weights -f w = Ylj^jj, 
we have more simplified expression for B a ^,~/(q,p) as 



B, 



cx.d 



l fom - 1 

,7(9. p) =^e e 

n=0 0^iuC/ s jGiu 
fc m -l s 

=-!+iEn 



1 + II i 1 + ^«,d(2/n,rf(j-l)+0) 



1 + 7jD ajd 

n=0 j = l 

Hence the result for the second part follows. 



-1 + ]J (l + <t> a ,d{Vn,d(j-l)+l)) 



1=1 



□ 
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4 Component-by-component construction of poly- 
nomial lattice rules 



4.1 Construction algorithm 

In order to efficiently find good polynomials q — (q±, . . . , qd s ) such that B at d.~f{q,p) 
becomes small, we investigate the CBC construction. We denote by Rb, m the 
set of all non-zero polynomials over ¥b with degree less than m, that is, 

Rb,m = {q € ¥b[x] : deg(q) < m and q ^ 0}. 

We search q from Rf s m component by component. Without loss of generality 
we can restrict q\ = 1. Thus, the CBC construction proceeds as follows. 

Algorithm 1 For a prime base b, a dimension s, an integer m, an interlacing 
factor d, and a smoothness a and weights 7 = ( r y u )u<ii a : 

1. Choose an irreducible polynomial p G F&[a;] such that deg(p) = m. 

2. Set qt = l. 

3. For r — 2,... ,ds, find q* by minimizing B ajC i^((q*_-i,q r ),p) as a function 
of q r £ Rb, m where q*_ 1 =(?£,..., q*_ x ). 

The following theorem gives a bound on B at d,~/(q*,p) for 1 < r < ds, which 
justifies the CBC construction. In the proof of the theorem, we shall use the 
following inequality that is sometimes referred to as Jensen's inequality. For a 
sequence (a n )nGN of non-negative real numbers 



where < A < 1. 

Theorem 2 Suppose q* ds = (ql, . . . , q^ s ) is found by Algorithm 1 for a prime 
base b, a dimension s 7 an integer m, an interlacing factor d, and a smoothness 
a and weights 7 = (7 tl ) u c/ a ■ Then for any r = 1, . . . , ds we have 



B aAy (q*,p) < 



(b m - 1)VA 



(6) 



1/A 



for l/(2min(a, d)) < A < 1, where we write r = d(jo — 1) + do such that 
jo, do € N and do £ {1, . . . , d} and 



G a ,d,\,a = D a d ( — 1 + (1 + Ca,d,A) a ) 




1 



' ^2Amin(a,d) 
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Proof. We prove the theorem by induction. For r = 1, we have q\ — 1 and 
Ba.,d,ii^,p) can be calculated as 

oo 

fe 1= l 

b m |fe 
oo 

= l{i}D a ,d (b - l)b l - m - 1 b- 2min ^ l - a 

l=m+l 

1 - 6-1 

— ^2min(Q,d)m^'{ 1 } Q ' d (^2 min(a,d) _ J,) 



< 



1 



(b m - 1) A 



7{l}G ! a,d,A, 



1/A 



Hence the result follows. Next we suppose that (6) has been established for 
r = d(jo — 1) + do. We then write r + 1 = d{j\ — 1) + d\ such that ji,d\ E N 
and d\ £ {1 , . . . , d} . It is obvious that 



(jo + 1, 1) if do = d, 
(jo, do + 1) otherwise. 



(h,di) -- 

We now consider 

B a ,d,y((q*,q r +i),p) 

= 7«(u)-Da,d E r Qjd (fc u ) 

(fc u ,o)eD- L ((q;, 9r+1 ), P ) 
= ^t.)^ 11 51 r a>d (k u ) 

(fc u ,0)G-D ± (q^p) 

+ T»(uU{r+l})^„,d {r+ X! ^,d(fc«,fc r+ i) 

(fc u ,fe r+ i)eNi"i +1 

(fc u ,A; r+ i,0)eD i ((q;,g r+ i),p) 

=B a4n {q* r ,p) + 6(q r +i), (7) 
where we have defined 

%r+i) := X! 7<uu{r+i})-D„,d + 1 5Z r a , d (k u , k r+ l). 

u ^ (fc„,Av+i)eN W+1 

(fc„,fc r+1 ,o)e_D- L ((q;, 9r+1 ), P ) 

In order to minimize i? Qj( j. T ((q*, as a function of qv+i, we only need to 
consider 9(q r+ i). Due to the averaging argument, the minimum value of Q(q r +i) 
has to be less than or equal to the average value of 6{q r+ i) over Rb. m - Applying 
Jensen's inequality, we have for l/(2min(a, <£j) < A < 1 



q r+ ieB, b/ 

<■ L_ ^ nA|«(«u{r+i})| 

~t,m_l Is 2s M«U{r+l})-^ a ,d 
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X 



(fc„,fe r+ l)(=Nl"l + 1 



(ku,k r+1 ,0)£D ± ((q' r ,q r+1 ),p) 

E X P,A|v(«U{r+l})| 
7-u(«u{r+l})- L 'a,d 



X 



E E ^,«l(*t„*r+l). 

<?r + iei?i,, m (fc u: fc r+ i)eNi"i+ 1 

(fc u ,fc r+1 ,o)e_D i ((q;, 9r+1 ), P ) 

For a fixed u C 7 r we have the following for the inner double sum in the last 
expression. If fc r+ i is a multiple of b m , we always have tr m (/c r +i) = and 
the corresponding term becomes independent of q r +i, or otherwise we have 
tr m (Av+i) ^ and tr m (fc r .+i)<7r+i cannot be a multiple of p by considering that 
p is irreducible. Hence we have 

-g^~l E E rl d (k u ,k r+1 ) 

q r+ l£R bim (fc u: fc r+ i)eNl"l + 1 

(k u ,k r +i,o)eD ± ((q*,q r+1 ), P ) 



oo 



= E r id(k r+ i) E r id(ku) 

6 m |fc r+ i tr m (fe„)-q u =0 mod p 

f) m ffc r+ i tr m (fc„)-q u ^0 mod p 

Here in the first term of the right-hand side we have 

Ea ( , \ b — 1 C a< d,x 
r a,d\ K r+l) — UXC2 min( a,d)m+a)(h2X mint a,d) _ h\ ~ 



£jA(2 min(a,d)m+a) ^2A min(a,d) ^ — ^2A min(a 

b""|fc r+ i 

and in the second term 

oo oo oo 

E r id(kr+l)= E r id(kr+l)- E r a,d(*V+l) 
A; r _j_l — 1 fcr+l — 1 &r-+i — 1 

6 m ffe r+ i 6 m |fc r+ i 

=li- 1 ^ ^ 



^2A min(a,d)m y ^aA ^2A min(a.d) ^ 

< C a ^d.,X- 

By inserting these inequalities into the above expression, we have 
E E r id(ku,K +1 ) 

q r+ ieR b , m (fc u ,fe,. + i)GN |u| + 1 

(k u ,k r +l,0)£D X ((q*,q T+1 ),p) 

^ Cq,rf,A \ A /l \ I Cg,d,X \ X /-jl \ 

-52Amin(a,d)m r a,dl K «J + ym _ ^ 2^1 r a,d[ K ^) 

fc„eN |u| fc u eN |u 

tr m {k u )-q u — mod p tr^ ) -q u ^0 mod p 
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n <7 H+1 



fc„GP 



Thus, we have a bound on x (q^i) as 

uCJ r 

We arrange the sum of the right-hand side of (8). We recall that r = (jo — 
l)d+do = (Ji— l)d+di — 1. Let us define two disjoint subsets J\ := {1, . . . , d(j\ — 
1)} and J2 := {d(ji — 1) + 1, ... , d(ji — 1) + d\ — 1}. In case of d\ = 1, J2 is 
empty. Every subset u C I r can be classified by whether it includes at least one 
element of J\ and by whether it includes at least one element of Ji- Since {r+1} 
is one of d components for the ji-th dimension, whether or not u includes some 
element of J 2 does not affect v(uU {r + 1}). Through this argument, we have 



<v(vU{r+l}) 1J a,d °a,<i,A 

iCI r 

E r \u 2 \ + 1 f ) \\v(u 1 )Ll{j 1 }\ r \u 1 \ fQ s 



U2QJ2 u 1 CJ 1 



We further arrange the inner sum of (9) by collecting the terms such that v(u\) 
u for u C For such terms, at least one component from {d(j — 1) 

1, . . . , dj} for every j £ u must be included. Thus we have 



E^AKuxMMl^l- 



P ) A|i;(«i)U{ii}| r ,|«i| 

/«(Ml)U{jl> 

«i C.7i 

= E 7 u u {jl} ^? l+1) [-l + (l + a )d ,A) d ] H 
For the outer sum of (9) we have 

E ^i^A 1 = Ca,d,\( l + C a ^\) dl ~ l 
u 2 CJ 2 

= (1 + C a ,d, x ) d i -(l + C a , d , x ) dl - 1 . 
By inserting these results into (9), we obtain 

E A f,A| t ,(«U{r+l})| r |«| + l 

«CJ r 

=(G a) d,A,di - Ga,d,X,di-l) X 7uU{ji } Gq ^ A (i 

From (8) we obtain a bound on x (q* +1 ). By using Jensen's inequality for (7), 
we finally have 



B a,d,-f( q *+fP) = (Bcd^iqUp) + 0(q* +1 ))' 
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1 



<- 



i 



T" G L"rf,A,d + G a,d,\,d Y ^uU{j } G a!d,\,d 

0#uC/ JO _! 



uCI 



30-1 



6™ - 

1 



— -J (C?a,d,A,di - G Qi d,A,di-l) ^ 7-uU{ji}G ! L,d,A,o 



6 m - 1 

Hence the result follows. 



"fu G a!d,\,d + G a,d,\,di 7«U{ji}Ga^ jA>d 



□ 



4.2 Fast component-by-component contruction 

Here we assume product weights for sake of simplicity and show how one can 
apply the fast CBC construction using the fast Fourier transform. The cost 
of the CBC construction by naive implementation of Algorithm 1 is at least of 
order 0(dsb 2m ) operations, which can be reduced to order 0(dsmb m ) operations 
for the fast CBC construction using the fast Fourier transform. 

According to Algorithm 1, we set q* = 1 and construct the polynomials 
Q21 ■ ■ ■ ilds inductively in the following way. Assume that q* = (q*,...,q*) 
are already found. Let r = d(jo — 1) + do and r + 1 = d(ji — 1) + d\ such 
that jo,do,ji,di G N and < do,d\ < d. As in the proof of Theorem 2, 
(ji,di) = (jo + 1,1) if d = d, or otherwise (ji,di) = (jo, do + 1)- Here we 
introduce the following notation 

d 

Vn,l ■= II I 1 7jA*,d + ljD a ,d J| (l + ^ a ,d(Vn,d(j-l)+l)) > 
j=l 1=1 



:= 

'In.r 



[ (l + <t>oi,d{yn,d{h-l)+l)) > 



and 

Vn.r : ~ Vn.rVn}-, 

(2) 

for < n < b m . We set rfk. r = 1 for < n < b m when d\ = 1, or equivalently 
when do = d. It is straightforward to confirm that we only need to minimize 

b m -l 

Vn,r<f>a,d(yn,r+l), (10) 

n=l 

for g r+ i G i?b /m . Since p is an irreducible polynomial over ¥b[x], there exists a 
primitive polynomial g in ¥{,[x]/p. When q r +\ = g % , (10) can be rewritten into 

b m -l 



£w^(^ (^)) 



18 



for 1 < i < b m . 

We now define the following circulant matrix 



t>a,d V 



Let us denote c = (ci, . 



p(a;) 

^Cbm-i) 1 " and »7 r = (r?i, r 

c = rt p r] r . 



l<i,n<b m 

i ■ • ■ , f?fe»"-i,r) T , and let 



Then for an integer io (1 < iq < b m ) such that Ci < for 1 < i < b m , we 
set q*. + i = g l ° . After finding we need to update rjn}- and rjn)- as follows. 

When di = d, 



„(i) 

'ln,r+l 



„(1) 



n {2) - 1 



1 A A (2) {tiff n ( x ) c ltj--i ( x ) 
l - J^D^d + 7j 1 I? a ,d?7n,r ( l + <P«,d K ( ^ 



Otherwise if di = l, . 




— >ln,r- 

_ ( 2 ) 

— i}n,r 



i 



»( a: )g*+l( :r ) 



Since the matrix f2 p is circulant, the matrix vector multiplication r2 p r/ T . can 
be efficiently done by using the fast Fourier transform as shown in [23, 24], 
which significantly reduces the computational cost as compared to the naive 
matrix vector multiplication. Whereas we focus only on product weights here, 
it is possible to apply the fast CBC construction to the case with another form 
of weights by minor modifications of the above procedure. We refer to [6, 17] 
for the fast CBC construction of lattice rules for order-dependent weights and 
POD (product and order-dependent) weights. 



4.3 Tractability properties 

Finally in this section, we briefly mention about the tractability properties of our 
algorithm. In the concept of tractability of multivariate integration, we study 
the dependence of the minimum number of points iV(e, s) on e and the dimension 
such that e {Q^u s \ s ,H s>(Xt ^) < ee 2 (Qo. s , -ffs.a.-y)- Given that the number of 
points is N — b m and that B a a t y(q,p) is a bound on e 2 (<3&™. s , H s-a ^), we have 
from Theorem 2 



N(e,s) < 70 



-X 



.d,X.d 



for 1/(2 min(a, d)) < A < 1. Hence, it is already obvious that N(e,s) depends 
polynomially on e" 1 with its degree at most one. As for the dependence on 
the dimension, we have the following corollary. Since the results immediately 
follow, we omit the proof. 

Corollary 2 Let b a prime, p G Vb[x] be an irreducible polynomial with deg(p) = 
m. Suppose that q = [qi, ... , q^s) is constructed according to Algorithm 1. Then 
we have the following 
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1. Assume that 



for some 1/(2 min(a, d)) < A < 1. Then N(e,s) is bounded above inde- 
pendently of the dimension. 

2. Assume that 



lim sup 



< OO, 



for some l/(2min(a, d)) < A < 1 and q > 0. Then N(e,s) depends 
polynomially on the dimension with its degree at most q. 



5 Numerical experiments 

Wc conclude this paper with numerical experiments. We present the values of 
B a d,-y{q, p) for different choices of s, a, d and 7, where q is found by Algorithm 1. 
In our computation, the base b is fixed at 2 and product weights arc considered. 
As a competitor, we use Sobol' sequences and Niederreiter-Xing sequences as 
interlaced components in Definition 5, instead of polynomial lattice rules. We 
use Sobol' sequences as constructed in [16] for any ds and Niederreiter-Xing 
sequences as implemented in [26] as long as 4 < ds < 16. For these rules, we 
denote a bound on the mean square worst-case error by 

B a ,d,-f(Ci, . . . ,Cds) = E 7v{u)D^d E r a> d{k u ), 

(k u ,a)eD- L (C u ...,c ds ) 

where C\ , . . . , Cds are generating matrices of the digital net and D (Ci , . . . , Cds ) 
is its dual space. 

We compare the values of B ai d y -y{q,p) with the values of -Ba^-^Ci, ■ • ■ , Cds) 
in Tables 1-4. In these tables, our constructed point sets based on polynomial 
lattice rules are denoted by PLR for short, and similarly, point sets based on 
Sobol' sequences and Niederreiter-Xing sequences are respectively denoted by 
Sobol' and N-X for short. In Tables 1 and 2, we consider 7^ = 1 for 1 < j < s, 
that is the so-called unweighted case. In Table 1, we fix a = d = 2 and change 
the dimension from s = 1 to s = 5. In Table 2, we fix s = 3 and change a and 
d simultaneously. In Tables 3 and 4, we do the similar comparisons for the case 
with 7j = j~ 2 for 1 < j ' < s. 

In most cases, PLR outperforms both Sobol' and N-X. A few exception can 
be seen for s = 3 and (a, d) = (1, 1) where Sobol' slightly outperforms PLR. We 
have to emphasize that what we compare here is not the mean square worst- 
case error itself, but a bound on the mean square worst-case error. Thus, it 
is possible that Sobol' and N-X outperform PLR in terms of the mean square 
worst-case error. However, our result is a good indicator of the usefulness of 
our algorithm in the sense that the actual performance of our constructed point 
sets is no worse than the result here. 
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Table 1: Values of B a ^,-y(q,p) and B a: d,-y(Ci, . . . , Cds) for 7^ = 1, 1 < j < s 



and (a, d) = 


[2, 2) with 


various 


choices of 


s = 1,2,5 








m 


s 


= l 


s = 2 


s = 5 




Sobol' 


PLR 


Sobol' 


N-X 


PLR 


Sobol' 


N-X 


PLR 


4 


2.13c-5 


2.11e-5 


2.69c-3 


2.56c-3 


2.70c-3 


1.27c+0 


1.76c+0 


9.81c-l 


5 


1.51e-6 


1.42e-6 


1.06e-3 


2.62e-4 


3.05G-4 


4.05e-l 


8.52e-l 


2.91c-l 


6 


1.38e-7 


9.56e-8 


9.51e-5 


2.04e-5 


7.58e-5 


1.84e-l 


1.50c-l 


7.42c-2 


7 


2.77e-8 


6.38e-9 


3.27e-6 


1.61c-6 


6.94e-6 


4.37e-2 


6.13e-2 


2.59e-2 


8 


l.llc-8 


4.24e-10 


4.34e-7 


1.74e-7 


4.82e-7 


2.15e-2 


1.84o-2 


6.55o-3 


9 


5.26e-9 


2.81C-11 


4.92e-8 


3.08e-8 


8.09e-8 


1.28e-2 


5.97o-3 


1.94o-3 


10 


2.62e-9 


1.86e-12 


1.32e-8 


l.llG-8 


5.78e-9 


9.43e-4 


4.54e-3 


3.97o-4 


11 


1.31c-9 


1.24e-13 


5.59e-9 


5.27c-9 


5.39e-10 


4.59e-4 


3.74c-3 


7.42c-5 


12 


6.55o-10 


6.44e-15 


2.57e-9 


2.61c-9 


4.64e-ll 


1.13e-4 


5.36e-5 


1.82c-5 


13 


3.28e-10 


4.44e-16 


1.28e-9 


1.30e-9 


4.85e-12 


8.07e-5 


1.05c-5 


4.32e-6 


14 


1.64o-10 


< 10~ 16 


6.37o-10 


6.52e-10 


3.99e-13 


1.38o-5 


1.62c-6 


7.18o-7 


15 


8.19C-11 


< 1CT 16 


3.19c-10 


3.19o-10 


4.35o-14 


7.52e-7 


2.39e-7 


1.35o-7 



Table 2: Values of B ay d :1 (q,p) and B ai d,-y{Ci, . . . , Cds) for jj = 1, 1 < j < s 



and s = 3 with various choices of (a, d) = (1, 1), (2, 2), (3, 3). 



m 


(a,d) 


= (LI) 


(a, d) = (2, 2) 


(a,d) = (3,3) 




Sobol' 


PLR 


Sobol' 


N-X 


PLR 




Sobol' 


N-X 


PLR 


4 


2.12o-2 


2.15o-2 


7.16o-2 


2.08e-l 


4.77o 


2 


1.86c+3 


1.79o+3 


1.14c+2 


5 


6.68e-3 


6.94e-3 


3.06e-2 


1.33c-2 


8.05e- 


3 


6.110+2 


7.22o+2 


1.87e+l 


6 


2.02e-3 


2.17o-3 


6.18e-3 


2.92e-3 


1.90e 


3 


2.58e+2 


4.06e+2 


1.14c+l 


7 


6.01e-4 


6.54e-4 


9.08e-4 


1.30o-3 


2.79e 


4 


2.07e+l 


1.01o+2 


1.35e+0 


8 


1.80e-4 


2.16e-4 


2.42e-4 


3.74o-4 


6.02e 


5 


3.55e+0 


9.83o+l 


1.34e-l 


9 


5.18o-5 


6.23e-5 


8.86e-6 


4.40e-6 


7.53e- 


6 


1.80e+0 


1.65o+0 


1.74e-2 


10 


1.50o-5 


1.86e-5 


1.58o-6 


7.80o-7 


9.00e 


7 


2.17c-l 


1.59o+0 


2.29e-3 


11 


4.29e-6 


5.50e-6 


1.20e-6 


1.66c-7 


1.45e 


7 


1.77e-2 


8.70e-3 


1.34o-4 


12 


1.21o-6 


1.61e-6 


6.41e-8 


1.69e-8 


1.61o 


8 


4.04e-3 


3.50e-3 


8.42e-6 


13 


3.39e-7 


4.43e-7 


7.57e-9 


4.92e-9 


3.08e 


9 


1.97e-3 


1.97e-3 


8.32o-7 


14 


9.53e-8 


1.34e-7 


2.43e-9 


2.06e-9 


2.37e- 


10 


9.82e-4 


1.13e-3 


5.14c-8 


15 


2.63o-8 


3.73o-8 


1.20c-9 


9.80o-10 


3.18o 


11 


4.91e-4 


4.79o-4 


2.75o-9 



Table 3: Values of B at d a (q,p) and B a ,d,-y(Ci, . . . , Cds) for jj = j 2 , 1 < j < s 



and (a, d) = 


(2, 2) with 


various 


choices of 


S = 1,2,1 








m 


s 


- 1 


s = 2 


s = 5 




Sobol' 


PLR 


Sobol' 


N-X 


PLR 


Sobol' 


N-X 


PLR 


4 


2.13o-5 


2.11c-5 


6.89c-4 


7.34c-4 


6.91e-4 


2.78c-2 


1.31e-l 


6.67e-3 


5 


1.51c-6 


1.42e-6 


2.66e-4 


6.76o-5 


7.72e-5 


3.33e-3 


1.04o-l 


1.38e-3 


6 


1.38e-7 


9.56e-8 


2.39e-5 


5.26e-6 


1.90e-5 


6.07e-4 


5.56o-4 


3.16e-4 


7 


2.77e-8 


6.38e-9 


8.38e-7 


4.30e-7 


1.74e-6 


1.61o-4 


1.32o-4 


6.41o-5 


8 


l.llc-8 


4.24e-10 


1.17c-7 


5.71c-8 


1.21c-7 


7.98e-5 


4.98e-5 


1.46o-5 


9 


5.26e-9 


2.81c-ll 


1.62e-8 


1.37o-8 


2.02e-8 


1.94e-5 


1.41o-5 


2.35e-6 


10 


2.62e-9 


1.86e-12 


5.27e-9 


4.74e-9 


1.45e-9 


2.27e-6 


3.53e-6 


5.09e-7 


11 


1.31c-9 


1.24e-13 


2.38e-9 


2.30e-9 


1.35e-10 


1.44e-6 


2.30e-6 


6.98e-8 


12 


6.55e-10 


6.44e-15 


1.13c-9 


1.14c-9 


1.16e-ll 


5.39e-8 


1.06o-7 


1.70o-8 


13 


3.28e-10 


4.44e-16 


5.65c-10 


6.91c-10 


1.21c-12 


3.38e-8 


1.63e-8 


2.69e-9 


14 


1.64c-10 


< lO" 16 


2.82c-10 


3.45e-10 


9.97e-14 


4.33c-9 


1.57o-9 


3.92e-10 


15 


8.19c-ll 


< lO" 16 


1.41c-10 


1.41c-10 


1.09c-14 


1.06c-9 


3. 91o-10 


7.29c-ll 
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Table 4: Values of B a ^^(q,p) and B a ,d,-y(Ci, . . . , Cds) for jj — j 2 , 1 < j < s 



and s = 3 with various choices of (a, d) = (1, 1), (2, 2), (3, 3). 



m 




= (1,1) 


(a, d) - (2, 2) 


(a,d) = (3,3) 




Sobol' 


PLR 


Sobol' 


N-X 


PLR 


Sobol' 


N-X 


PLR 


4 


3. 74c- 3 


3.67c-3 


5.49c-3 


2.21c-2 


2.38o-3 


1.12o+2 


1.310+2 


6.13O+0 


5 


1.04e-3 


1.04e-3 


2.03e-3 


9.73c-4 


4.25e-4 


1.86e+l 


2.42o+l 


6.03e-l 


6 


2.85e-4 


2.85e-4 


2.00e-4 


2.84o-4 


9.00e-5 


7.21e+0 


1.13o+l 


3.720-1 


7 


7.77o-5 


8.04e-5 


2.84e-5 


3.93e-5 


1.37e-5 


5.82e-l 


2.81o+0 


5.32e-2 


8 


2.16e-5 


2.23e-5 


6.95e-6 


1.14c-5 


2.21e-6 


1.03c-l 


2.73o+0 


4.580-3 


9 


5.74e-6 


6.04e-6 


2.87e-7 


2.64o-7 


2.53o-7 


5.10e-2 


4.60e-2 


5.02e-4 


10 


1.58e-6 


1.68e-6 


5.04e-8 


7.85e-8 


3.22e-8 


6.09e-3 


4.41c-2 


7.550-5 


11 


4.25e-7 


4.71c-7 


3.61c-8 


8.29e-9 


4.35e-9 


5.20e-4 


2.67e-4 


3.98e-6 


12 


1.14e-7 


1.27e-7 


3.01c-9 


1.93e-9 


5.93e-10 


1.26o-4 


1.10e-4 


2.38e-7 


13 


3.05e-8 


3.37e-8 


8.19G-10 


8.93e-10 


9.78e-ll 


6.17o-5 


6.10o-5 


2.32e-8 


14 


8.25e-9 


9.52e-9 


3.72o-10 


3.79e-10 


7.46e-12 


3.08e-5 


3.39e-5 


2.02e-9 


15 


2.17o-9 


2.52o-9 


1.85e-10 


2.27o-10 


1.14o-12 


1.54o-5 


1.520-5 


1.05o-10 
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